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We report first thermodynamic measurements of the temperature derivative of chemical potential {dfi/dT) 
in two-dimensional (2D) electron systems. In order to test the technique we have chosen Schottky gated 
GaAs/AlGaAs heterojunctions and detected experimentally in this 2D system quantum magnetooscillations 
of dp/dT. We also present a Lifshits-Kosevitch type theory for the dp/dT magnetooscillations in 2D systems 
and compare the theory with experimental data. The magnetic field dependence of the dp/dT value appears to 
be sensitive to the density of states shape of Landau levels. The data in low magnetic field domain demonstrate 
brilliant agreement with theory for non-interacting Fermi gas with Lorentzian Landau level shape. 


Quantum oscillations are known to be a universal tool 
to study the electron energy spectrum (Fermi surface 
cross-sections, electron effective mass and g-factor) in 
three-dimensional single crystals and two-dimensional 
(2D) systems. In contrast to the 3D-case, the 2D sys¬ 
tems allow in-situ tuning the spectrum and the Fermi 
energy by various methods (including electric field ef¬ 
fect in gated structures, illumination, uniaxial stress 
etc.), and, hence, allow comprehensive magnetooscilla¬ 
tion studies. Quantum oscillations in 2D systems are 
most often studied in resistivity (Shubnikov-de Haas ef¬ 
fect) [T], and magnetization (de Haas-Van Alphen) [5]. 
However, other physical quantities in 2D systems, such 
as thermo-EMF [3], heat capacity [2], chemical poten¬ 
tial EE], compressibility [7], also oscillate with mag¬ 
netic field. All experimental methods applied for mea¬ 
suring magnetooscillations of various observables have 
limitations for a practical use, therefore developing new 
magnetooscillations tools is of substantial importance. 

All methods for chemical potential measuring are 
based on the W. Thomson (Kelvin) idea [5] that if the 
two plates of a capacitor are made of different materials 
(with different work functions Ri, 2 ), the charge of the 
plates is C{ri — where C is the electric capaci¬ 

tance. Correspondingly, when the two plates are con¬ 
nected electrically and an external parameter varies af¬ 
fecting one of the chemical potential values, a recharging 
current starts flowing between the plates. The recharg¬ 
ing current is proportional to A/i (in case the capaci¬ 
tance C varies), oc dfi/dn (in case one of the plates is a 
2D gas of a density n varying with a gate voltage m\), 
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OC dfL/dB (in case the magnetic field B varies [U]), or 
oc dft/dT (in case the temperature varies [TU]). 

In our study we apply the technique of measuring 
dfi/dT similar to that used by Nizhankovskii for bulk 
samples m, to the 2D electron systems in magnetic 
field. We focus on dfi/dT oscillations to compare them 
with semiclassical theory, and to determine the shape of 
the density of states of Landau levels. The advantage of 
the temperature modulation technique is the absence of 
eddy currents or a background signal (concomitant of 
many other techniques), and its pure thermodynamic 
origin. 

Qualitative discussion. - It is worthwhile to give a 
qualitative explanation why dft/dT oscillates with per¬ 
pendicular magnetic field. For the bare quadratic en¬ 
ergy spectrum, e{p) = p‘^/2m, the single electron den¬ 
sity of states is constant in two dimensions. At temper¬ 
atures T Ep the number of particle-like excitations 
above p equals to the number of hole-like excitations be¬ 
low p (hatched areas in Fig. [T^). Therefore, for a fixed 
electron density n the chemical potential is indepen¬ 
dent of temperature, dp/dT — 0, with exponential ac¬ 
curacy at low temperatures, T <C Ep, where Ep stands 
for the Fermi energy. In the case of energy dependent 
density of states (like e.g. in 3D systems, graphene or 
2D systems in quantizing magnetic field), one can ex¬ 
pand it in the vicinity of the Fermi energy (see FiglUD): 
D{e) = D{Ep) + {dD{e = Ep)/de) x {e - Ep) Then, 
for a nonzrero temperature the particle-hole asymme¬ 
try emerges (since the hatched areas do not coincide in 
FigHb) and one needs to shift chemical potential in or¬ 
der to conserve the total number of particles. Using the 
standard low temperature expansion for the Fermi-type 
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(a) Density of states D(e) = const and its prod¬ 
uct with the Fermi distribution function. Equivalent 
hatched areas denote electron an hole excitations, (b) 
The same as (a) for the density of states changing 
with energy; the hatched areas are unequal, (c) Os¬ 
cillatory density of states in perpendicular magnetic 
field and Fermi function corresponding to low temper¬ 
ature (T < hu!c,r, solid line) and high temperature 
(T < hujc,!', dashed line). 


integrals (see, e.g., Ref. [TT]) one can easily find correc¬ 
tion to the chemical potential at low temperatures, T 
Er: fi(T) - Ep = /6 D{Ef)] x dD{e = EF)/de. 

Hence, we find 

tt^T dDiEp) -K^TdD{EF) 

We note that this equation is applicable for degenerate 
Fermi-systems of any dimensionality and for any spec¬ 
trum. 

In perpendicular field, due to Landau quantization 
D{e) becomes dependent on energy (see Fig. [TJ;), and 
Eq. O should be averaged over T in the vicinity of the 
Fermi energy. If the temperature is low (Ti, Fig. [TJ:), 
one can directly apply Eq. (El); in the opposite limit of 
high temperature(r 2 . Fig. [TJ:) the oscillations are aver¬ 
aged over a wide energy interval and become exponen¬ 
tially suppressed. 

Two qualitative predictions follow from the above 
considerations. (i) In the maxima and minima of 
the density of states the dti/dT signal is zero, and 
vice versa, the signal is maximal where the derivative 
d\r\D{EF)/dEF is maximal, (ii) The amplitude of the 
dp-/dT magnetooscillations is a non-monotonic function 
of temperature: at low temperatures T <C WcF the 
dp/dT oc T, whereas for high temperatures T > lOc 
averaging over several oscillations suppresses the signal. 

Theoretical background. - The thermodynamic po¬ 
tential of non-interacting electron system can be written 
as [TT] 


n { T , p , B ) = -T 


deD{e) In 


I -b 


( 2 ) 


In the presence of the perpendicular magnetic field B 
the density of states becomes 

OO 

(T=± n=0 


where Z = gLpBB/2 {gp stands for the g-factor) de¬ 
scribes the effect of the Zeeman splitting. The func¬ 
tion >V(e) describes broadening of a Landau level due 
to disorder. It satisfies the normalization condition: 
/ deyV{e) = I. We note that, in general, this func¬ 
tion can be different for different Landau levels. Using 
the function W(e) a typical width of a Landau level can 
be estimated as F ~ [J de>V"(e)] . The accurate 

quantum-mechanical evaluation of yV(e) for a given type 
of disorder is a complicated problem and solved only 
partially [HHUj. In the absence of disorder, one obvi¬ 
ously has yV{s) = 6{e). There are three commonly used 
models of the disorder Landau level broadening [HUS]: 


Lorentzian model: 
Gaussian model: 
semicircle model: 


>V(£) = 


1 F 

TT £2 -b F^ 


>V(£) 



>V(£) 
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In a standard way, by means of the Poisson resum¬ 
mation formula applied to the thermodynamic potential 
(ED we obtain the Lifshitz-Kosevich-type [20] expression 
for dp/dT: 
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This result holds for ^ ± Z ^ WcT, F. Here X^ = 
27r^Tfc/a;c, and 

Ak= J deWie) exp (6) 

characterizes damping of oscillations due to Landau 
level broadening. For three models in Eq. (jT] it be¬ 
comes 


Lorentzian model: 
Gaussian model: 

semicircle model: 


Ak = 

7^2 1 2 / 2 
Ak = e-" ^ 


Ak 


UJc j f 27:Tk 
'KVk ^ \ OJc 


( 7 ) 
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where Ji(x) denotes the Bessel function. 

For low temperatures, T <C cUc, B, /x ± Z, the 
dp/dT value is given by the result ([T]) with D = 
(1/2) X]ct=± Hsre Ep = 2'Knr,lm denotes the 
Fermi energy for a a given spin projection. The corre¬ 
sponding electron density is determined by the Zeeman 
splitting, ria jn = 1-1- aZjEp where Ep = 2'Knlm. In 
the case of the Lorentzian broadening of a Landau level 
the density of states at the Fermi energy becomes 

D{Ei) = " (l/2)snl.(2„r/a,,)- 

TT sinh^(7rr/a;c)-b cos2(7r£l^/a;c) 

We mention that above we ignore the effect of 
electron-electron interaction. As well-known, for an 
interacting 3D electron system the Lifshitz-Kosevich- 
type expressions for magnetooscillations of the ther¬ 
modynamic potential can be obtained from the non¬ 
interacting expressions via standard Fermi-liquid renor¬ 
malization of the quasiparticle spectrum Eillig. In 
2D, it is not the case, in general [53]. For classically 
weak perpendicular magnetic held, modihcation of the 
Lifshitz-Kosevich-type expression for a 2D electron in¬ 
teracting disordered system has been studied in Refs. 

1531155]. 

Experimental results. - We studied two gated 
GaAs/AlGaAs single hetero junctions (GaAsl and 
GaAs2, similar to tat from Ref [^) in the temperature 
range between 2.5 and 25K with electron mobilities 20 
and 25 m^/Rs, respectively. The modulation technique 
was similar to that used by Nizhankovskii m and will 
be reported elsewhere [57]. The density could be varied 
in the interval from 2 • I0^^cm“^ to 4 • I0^^cm“^ by 
changing the gate voltage. The samples had 5mm^ 
area, the gate-to-2D gas capacitance was 1100 pF. 
Both samples demonstrated similar results, we present 
therefore the data for GaAs2 only. 

Examples of low-held dp/dT oscillations measured for 
the sample GaAs-2 are shown in Fig. [S] For helds lower 
than 3 Tesla and in the range of temperatures 2.7-9.IK 
we htted the data using Eq. (©. We use the bare band 
mass {mb = 0.067me for GaAs). Because of the low- 
held range studied and large cyclotron splitting, all the 
results below are insensitive to the g-factor value (that 
may vary in the range g = -0.4-^ 2). The best ht is ob¬ 
tained for the broadening of Landau levels described by 
the Lorentzian model with the width F=0.4meV which 
is independent of B (see dashed curves in Fig. 131). The 
range of temperatures (0.22mV< T <0.8 mV) and mag¬ 
netic helds (0.86 < hwc < 5.2 mV) used for htting pro¬ 
cedure is wide enough to identify the line shape with 
the Loretzian curve. On the other hand, the Gaus¬ 
sian model with F independent of B (solid curves in 
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Fig. [2] dp/dT versus magnetic held in sample GaAs- 
2 at various temperatures, n = 4.16 • 
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Flg-O dp/dT versus magnetic held in sample GaAs- 
2 at n = 4.16 • 10^^cm“^. Data is htted with 
Lorentzian (dashed curves),Gaussian (solid curves) 
and Gaussian with F oc '/B (dotted curves) line- 
shapes at 3K (panel a) and 9K (panel b). Upper 
axis shows hlling factor. 

Fig. [3]) fails to ht the oscillations. The Gaussian model 
with level broadening proportional to \/B works only 
in low magnetic helds (dotted curves in Fig. |31). We 
checked that in the density range from 3T0^^cm“^ to 
5T0^^cm“^ the low-held level broadening is constant 
(equal to 0.4 meV for sample GaAs-2). The Lorentzian 
lineshape of the Landau levels is in the agreement with 
magnetization measurements of Potts et al. m on the 
moderate mobility GaAs-based sample. Generally, the 
shape of DOS in 2D systems is still debated, with a 
history of research reviewed, e.g., by Usher et al. [5]. 
Different approaches are used to determine DOS with 
many contradictory results obtained. 

Conclusions. - We presented here experimental test 
of the novel technique of thermodynamic dp/dT mea¬ 
surements, that is ideally suitable for 2D gated carrier 
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systems. The technique appears to be highly sensitive 
and thereby enables thermodynamic measurements with 
a single-layer electron system comprising only 10^° elec¬ 
trons. We also present a Lifshits-Kosevitch type calcula¬ 
tions for the dfi/dT magnetooscillations in 2D systems 
and compare the theory with experimental data; the 
comparison reveals a good agreement between the data 
and the theory. The magnetic field dependence of the 
dp,/dT appears to be rather sensitive to the shape of the 
density of states at Landau levels, and in the particular 
case of moderate mobility samples studied we found the 
density of states to correspond to the Lorentzian rather 
than Gaussian curve. We want to emphasize, that the 
method described here suggests an independent oppor¬ 
tunity of extracting energy spectrum of the 2D carrier 
system. 
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